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Abstract 

Gallium chalcogenides are promising building blocks for novel van der Waals heterostruc¬ 
tures. We report on low-temperature micro-photoluminescence (PL) of GaTe and GaSe films 
with thicknesses ranging from 200 nm to a single unit cell. In both materials, PL shows dra¬ 
matic decrease by 10 4 -10 5 when film thickness is reduced from 200 to 10 nm. Based on 
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evidence from cw and time-resolved PL, we propose a model explaining the PL decrease as 
a result of non-radiative carrier escape via surface states. Our results emphasize the need for 
special passivation of two-dimensional films for optoelectronic applications. 


Introduction 

The discovery and research of remarkable properties of two-dimensional (2D) sheets of carbon, - 
known as graphene, has sparked interest in other layered materials such as metal chalcogenides 
(MCs). 2-3 Recent progress in device fabrication opens up new possibilities for 2D MC films in 
nano- and opto-electronics.- s Recent efforts in research of 2D semiconductors and their devices 
focused mainly on molybdenum and tungsten dicha 1 cogcnides--~ 13 and relatively little has been 
reported on optical properties of thin films of III-VI materials: only recently electronic- 1 ^ -16 and 
photonic^ trial devices have been reported for GaSe and GaTe, 16 and room-temperature photo¬ 
luminescence (PL) was measured in detail for InSe thin films. 18 As the crystal growth and im¬ 
provement of thin GaSe films continues, 14 ’ 19 ’ 20 further exploration of physical properties of III-VI 
materials presented in this work is motivated by their potential use as building blocks for novel 
van der Waals heterostructures, where materials with a range of bandgaps and band offsets will 
be required.-^ Furthermore, in contrast to molybdenum and tungsten chalcogenides emitting light 
efficiently only in films with a single unit cell thickness, III-VI materials are bright light emitters in 
a range of thicknesses. 18 This may relax the stringent fabrication requirements and add flexibility 
for the novel heterostructured devices such as light emitting diodes. 8 

Both GaTe and GaSe are layered crystals with strong covalent in-plane inter-atomic bonding 
(with some ionic contributions 21-22 ) and weaker predominantly van der Waals inter-plane bond¬ 
ing. 21-2326 Single tetralayer having hexagonal in-plane structure consists of two Ga atoms and 
two Se or Te atoms: Se-Ga-Ga-Se and Te-Ga-Ga-Te. 25 The bulk lattices are built by stacking 
tetralayers, which can occur in several ways. 21-23-24 For the wider studied GaSe, several types 
of stacking exist leading to different polytypes.' 21-23-24 For GaTe having a monoclinic crystal lat- 
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tice 22,27 the polytypic behavior has not been observed. 22,26 This may lead to lower probability of 
stacking faults in GaTe resulting in clearer observation (compared to GaSe) of excitonic features 
in optical experiments. 22,26 

Here we study optical properties of GaSe and GaTe thin films as a function of the film thickness. 
Continuous-wave (cw) and time-resolved low-temperature micro-PL for a wide range of films from 
200 nm to one tetralayer thicknesses is measured. PL intensity is used to monitor the quantum yield 
(QY), which falls dramatically for thin films: integrated cw PL intensity drops by up to ~ 10 4 — 10 5 
when the film thickness is reduced from 200 to 10 nm. A similar observation of reduced PL for thin 
films was previously reported for InSe and was explained as transition to a band-structure with an 
indirect bandgap, the conclusion also based on the observed PL blue-shift with the decreasing film 
thickness.— Such indirect bandgap behavior is also theoretically predicted for single monolayers 
of GaSe and GaTe. 25 However, no size-quantization effects as in InSe are observed in our work for 
GaSe and GaTe. Based on the evidence from both cw and time-resolved spectroscopy, we develop 
a model that shows that the observed PL reduction can be explained by non-radiative carrier escape 
to surface states. Our explanation does not require introduction of the direct-to-indirect band-gap 
transition. Following the cw PL data analysis, we identify a critical film thickness of about 30-40 
nm, below which the non-radiative carrier escape changes its character. The importance of surface 
states predicted by our results emphasizes the need for development of novel surface passivation 
for III-VI films, possibly involving oxygen-free dielectrics such as boron nitride. 28 

Experimental procedure 

Fabrication of GaTe and GaSe samples 

Single crystals of GaSe and GaTe were grown by high-pressure vertical zone melting in graphite 
crucibles under Ar pressure. The detailed description of crystal growth processes and properties of 
GaSe and GaTe can be found in Refs.—~— The gallium mono-chalcogenides used in this work were 
synthesized from high-purity materials: Ga and Te - 99.9999 %, Se - 99.9995 %. The crystals used 
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are high-resistivity semiconductors with low free carrier absorption, which has been confirmed in 
infra-red transmition measurements. GaSe and GaTe have n-type and p-type conductivity, respec¬ 
tively. This is a typical observation: selenides are usually of n-type conductivity, whereas tellurides 
often can have conductivity of both types, even within one ingot. Such behavior is attributed to 
deviations of the crystal composition from stoichiometry, usual for metal chalcogenides. In our 
case the n-type conductivity clearly indicates some excess of Ga (donor) in GaSe, whereas the 
p-type conductivity indicates a slight excess of Te (acceptor) in the GaTe. 

The III-VI thin films studied in this work were fabricated by mechanical cleaving from bulk. 
The films were deposited on Si/Si02 substrates. Within the first 15 minutes after the exfolia¬ 
tion/deposition procedure, the films were placed in a plasma-enhanced chemical vapor deposition 
(PECVD) reactor and a 15 nm S^INU layer was deposited with the sample maintained at a tem¬ 
perature of 300 °C. This process leads to a complete encapsulation of the films, protecting them 
from interaction with oxygen and water present in the atmosphere. Films with a wide variety of 
thicknesses were obtained from single unit cell (single monolayer, ML) shown in Fig jT|to 200 nm. 
The thicknesses of the films were determined using atomic force microscopy. 

Optical characterization methods 

Optical characterization of the GaTe and GaSe thin films was carried out using a low-temperature 
micro-photoluminescence (/iPL) technique. The sample was placed on a cold finger in a continu¬ 
ous flow He cryostat at a temperature of 10 K. A microscope objective with NA of 0.6 was placed 
outside the cryostat and was used to focus the laser beam on the sample (in a ^2/dm spot) and to 
collect the photoluminescence (PL) from the films. In continuous-wave (cw) experiments, PL was 
detected with a 0.5 m spectrometer and a liquid nitrogen cooled charge coupled device. For cw 
PL excitation a laser emitting at 532 nm (2.33 eV) was used. In the ultra-fast spectroscopy exper¬ 
iments the excitation of GaTe and GaSe layers was performed using frequency-doubled titanium- 
sapphire (wavelength of 415 nm, pulse duration ^ 200 ps) focused on the sample in ?»10/im spot. 
Time-resolved PL (TRPL) was detected using a streak camera. The temporal resolution of the 
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experimental setup for the time-resolved measurements was 10 ps. 


Experimental results 

Low-temperature cw PL results: GaTe 

Figj2[a) shows typical PL spectra measured for S^Nzi-capped GaTe thin films at T — 10K (cw 
laser power P =2 mW in FigJIJ'a)). In this figure, in all films but the one with the thickness hfu m =8 
nm, a narrow feature is observed around 1.75 eV. It is observed at an energy where free exciton 
(FE) PL is expected. We will therefore refer to such features in GaTe (and GaSe) films as a ’free 
exciton’ peak as opposed to the low energy broad PL bands corresponding to excitons bound to 
impurities/defects and observed in the range of 1.6-1.7 eV for GaTe. The ’free excitons’ may also 
experience disordered potential and degree of localization as evidenced from a relatively broad FE 
PL line of 10-15 meV (varying from sample-to-sample). A different behavior of the FE peak com¬ 
pared to the bound excitons has been verified in temperature and power-dependent measurements 
and is further confirmed in time-resolved studies discussed below. The FE peak is usually more 
pronounced in thick films of around 100 nm and above. In thin films as in the 8 nm film in the 
figure, the FE peak could only be observed under high power pulsed excitation, when the impu¬ 
rity/defect states saturate. The sharpest PL feature in GaTe spectra, the FE line, has the peak energy 
varying from film to film in the range 1.74-1.76 eV showing insignificance of size-quantization ef¬ 
fects in contrast to InSe where the PL blue-shift was observed for thin layers. 18 

A pronounced feature of PL measured from different films is a dramatic decrease of PL in¬ 
tensity with the decreasing thickness of the material [see Figj2[b)]: about 10 5 (10 4 ) decrease 
is observed between 200 (100) and 7 nm. The strongest PL reduction by 3 orders of magni¬ 
tude is detected between 200 and 40 nm. For hfu m <40 nm, the PL intensity reduction slows 
down and decreases less than 100 times when hjn m is varied between 40 and 7 nm. The dot¬ 
ted curve in the graph shows the expected PL intensity behavior assuming thickness-independent 
quantum efficiency, i.e. when reduction in PL is caused solely by the reduced absorption and 
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reduced number of e-h pairs created by the laser. The curve is described by the expression 
Ipl = lGaTe[\ ~ exp(-a GaTe hfn m )], where l Gate is the PL intensity for films with h f u m -200 nm. 
The absorption coefficient CCG a Te =5000 cm -1 is used according to Ref.- A discrepancy by a few 
orders of magnitude between the experiment and the calculated curve in a wide range of film 
thicknesses is evident on the graph. 

Low-temperature cw PL results: GaSe 

FigJUa) shows typical PL spectra measured for Si 3 N 4 -capped films of GaSe of several thicknesses 
between 8 and 70 nm (the cw laser power of 2 mW is used). The PL signal is observed in a 
range from 1.95 to 2.05 eY, which is below the emission energy of the free exciton, reported to 
be around 2.10 eV for some high purity GaSe samples (see e.g. Ref. 34 ). The detected PL in our 
samples thus comes from impurity/defect states. The observed localized states may originate from 
the non-stoichiometric composition of the bulk material. 

In Fig{3ja) it is observed that PL spectra of thin films <20 nm usually consist of multiple 
pronounced lines [a feature similar to GaTe in Figj2}a)], whereas PL spectra tend to exhibit a 
single pronounced peak for thicker films. PL linewidths vary between 15 and 60 meV. In some 
GaSe films, sharp PL lines with linewidths below 5 meV, similar to the sharp features observed in 
the PL spectrum for the 8 nm film in Figj3^a), also occur in the whole energy range of GaSe PL. 
From the data measured on more than 50 films, we observe that PL peak energies have a very wide 
distribution in the range 1.99-2.06 eV. Similarly to GaTe, there is no evidence for size-quantization 
effects as a function of the film thickness. 

Further evidence for exciton localization in thin films is a pronounced Stokes blue-shift ob¬ 
served when the laser excitation density is increased as shown in Fig 01 The inset in Figj4ja) 
shows that the PL peak shifts by ~20 meV as the cw laser power is changed from 0.01 to 2 mW. 
This is a typical behavior observed in all GaSe films independent on the film thickness: at high 
power, saturation of some of the PL features is observed accompanied in most cases with a blue- 
shift of PL of around 10-20 meV. This is a typical behavior observed for localized exciton states in 
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semiconductors, the effect also similar to the state-filling phenomenon in semiconductor quantum 
dots. 35 In some GaSe films, if the optical pumping is further increased, for example, by using 
pulsed excitation, a relatively broad free exciton feature can be observed, as shown in FigBlb). 
This behavior is in agreement with that observed previously in GaTe and GaSe under pulsed ex¬ 
citation, and is related to saturation of the localized states having relatively slow recombination 
rates. 36 

Similarly to GaTe films, a significant decrease of PL intensity with the decreasing thickness 
of the GaSe films is observed [see Fig[3tb)] by about 2x 10 4 between 200 and 7 nm. As for the 
GaTe films in Figj2[b), the strongest PL reduction by 3 orders of magnitude is detected between 
200 and 30 nm. For hfu m <30 nm, the PL intensity reduction slows down and is about 30 when 
hfu m is varied between 30 and 7 nm. Similar to GaTe in films with hju m <7 nm PL is completely 
suppressed. Similarly to Figj2] we show a curve that describes PL reduction due to the reduced 
absorption only calculated as Ipp — IcaSe [ 1 — ex P (— (%GaSehfilm)] , where IcaSe is the PL intensity for 
films with h fn m = 150 nm and the absorption coefficient a= 1000 cm -1 .- As for GaTe, a significant 
discrepancy by a few orders of magnitude between the experimental results and the calculated 
curve is clear in a wide range of film thicknesses. 

Time-resolved PL measurements 

In order to shed further light on the results of cw PL and also to provide further experimental 
foundation for our theoretical model, time-resolved PL experiments have been carried out for a 
range of films with different thicknesses. Fig[5] shows typical TRPL data obtained at T ~10K. 

For GaTe, the difference in the origin of the PL features observed in Figj2[a) is further evi¬ 
denced in TRPL. Fig[5]>i.b) shows data for a 27 nm thick film exhibiting a behavior typical for 
films with hfii m in the range 20 to 200 nm. Figj5[a) presents a streak-camera scan clearly showing 
two pronounced features at 1.76 eV and 1.71 eV corresponding to the free and localized excitons, 
respectively. The free exciton peak intensity decays considerably faster than that of the localized 
states, which does not change significantly on the time-scale of 130 ps shown in the figure. FigJ51b) 
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shows two decay curves measured at 1.76 eV (gray) and 1.71 eV (red). Fitting with single expo¬ 
nential decay functions is shown with blue curves and gives 10 ps for the free exciton and 150 
ps for the localized states. The inset shows a PL decay curve measured at 1.71 eV on a larger 
time-scale, exhibiting an almost complete decay of the signal at 1 ns. We find similar life-times 
for other films, however no clear dependence on the film thickness is observed: the free exciton 
PL decay time varies between 5 and 25 ps and that for the localized states between 100 and 200 
ps. The lifetimes also weakly depend on the laser excitation power. We also note the rise times of 
«H5 and 20 ps for the free and localized states, respectively, indicating fast carrier relaxation into 
the light-emitting states. 

Similar difference between the PL dynamics of the high and low energy part of the spectrum is 
observed for GaSe thin films in Figj5^c),(d). Here a behavior resembling the Stokes shift shown in 
Figj4] is observed: as the carrier density decreases with time after the laser pulse, the PL intensity 
maximum progressively moves to lower energy. Figj5jd) details the behavior shown in FigJ5}c): 
two decay curves measured at 2.055 eV (gray) and 2.035 eV (red) are shown. In the center of the 
PL band at 2.035 eV, the non-exponential decay occurs with a characteristic time of 400 ps, which 
also shows a slow-decaying component. At around 2.055 eV, the initial PL time-dependence can 
be well fitted with a single-exponential decay with a lifetime of ~40 ps, dominated most likely by 
carrier relaxation to lower energy. The complex behavior in GaSe films occurs due to the partial 
saturation of the states at short times after the excitation pulse and fast relaxation to lower energy. 
We find similar behavior for films with other thicknesses. Similarly to GaTe no clear dependence 
on the film thickness is observed and the rise times of f»30 ps are found for the localized states. 

Modeling 

In order to describe the observed trend of PL intensity as a function of film thickness hfu m in GaSe 
and GaTe thin films we have developed a rate equation model presented in detail in the Supplemen¬ 
tary Information (SI). 37 As shown in Figjbja), in the model we consider the following processes 
leading to light emission: optical excitation of e-h pairs (uniformly across the full thickness of the 
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film); relaxation into the non-radiative traps with a time T nr \ or z nr 2 depending on the thickness 
of the film, as explained below; relaxation with a time T,- e / into the light-emitting states denoted 
in Figjfja) as ’PL states’; PL emission from these states with a time Tpp. The absence of a clear 
dependence of the PL decay times on the film thickness allows us to use a thickness-independent 
constant T pp and neglect any other decay mechanisms for the e-h population of the light-emitting 
states. The physical origin of this may be in a relatively strong localization of e-h pairs in the 
light-emitting ’PL states’, which leads to suppression of non-radiative processes. We also assume 
that the mechanism of relaxation into the ’PL states’ is not sensitive to the film thickness and can 
be described with a thickness independent time T,- e /. 

As shown in Figj6}b) and (c) we also assume that the film is divided in three regions: (1) 
two regions of thickness ho near the film surfaces where fast carrier relaxation to surface traps 
occurs leading to non-radiative carrier escape; (2) a ’normal’ region of thickness hfu m — 2/?o in 
the film’s central part not containing the traps where the photo-excited carriers can escape into 
regions (1), where they undergo non-radiative decay. Here the ho value may be associated with a 
depletion depth or an average surface trap radius.- As shown in Figjb] there are two possibilities: 
in FigHtb), where hfu m > 2 /zq and both regions of type (1) and (2) exist; and in Figjfjc), for films 
with hfu m < 2/i(), where region (2) is not present. The model assumes that light absorption and 
PL occurs in both types of regions. In both regions, e-h pairs relax with the time T re / into the 
light-emitting states. 

In region (1), in first approximation the average time it takes for the carrier/e-h pair to escape 
non-radiatively is proportional to half the thickness of region (1) (can be understood as the average 
time for the carrier to reach the surface or as the overlap of the wavefunction of the carrier and the 
surface trap). In region (2), the non-radiative escape time reflects the average time it takes for a 
carrier or an e-h pair to reach any of the regions (1). The underlying mechanism for this process 
may be depletion and band-bending expected at the film surface leading to charge separation and 
non-radiative decay. 1 We assume that once the carrier or e-h pair has reached region (1) it escapes 
non-radiatively. In first approximation, the average time it takes a carrier/e-h pair to reach region 
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(1) is proportional to half the thickness of region (2). Thus we introduce non-radiative decay times 
in region (1) as x nrl = (h 0 /2)/u x for h fUm > 2 h 0 , and x nr \ = (hfu m /4)/u\ for h fibn < 2 h Q . In 
region (2) it is x nr 2 — {hfu m /2 — ho)/u 2 ■ Here u x and 112 have dimensions of m/s. In the case of 
region 2 where effectively we assume ballistic exciton (or electron/hole) transport preceding the 
non-radiative escape, 112 can be interpreted as the average carrier velocity. 

Discussion 

As detailed in Supplementary Information 37 and observed in Figj2|b) and Fig{3}b), we find that the 
proposed model provides a reasonable description of our data using four parameters (which are not 
completely independent as we find): the ’critical’ thickness ho, a parameter describing the amount 
of light absorbed by the film, and products x re iu\ and x re /U 2 . In particular the fitting functions that 
we produce capture the change in the ’slope’ of the data observed at around 30 nm for GaSe and 40 
nm for GaTe, which we interpret as the thickness of the film where hju m & 2/?o, i.e. the thickness 
of region (2) turns to zero, and non-radiative carrier escape changes its character. 

However, we find that the accuracy of the fitting is not sufficient to extend our analysis beyond 
determination of the order of magnitude of the products x re iu\ and x re iU 2 - We find that x re iu\ 
and T,- e /M 2 are of the order of 1000 nm and 100 nm, respectively, for both GaSe and GaTe. The 
solid line shown in Figj2](b) for GaTe films is obtained for x re [U X = 3800 nm and x re ]U 2 — 200 nm, 
whereas the fitting in Figj3](b) for GaSe films is done for x re iu x = 2900 nm and x re [U 2 = 150 nm. 
We note that the description of GaSe PL is more satisfactory, possibly because in GaTe there is 
a contribution from free exciton PL, so additional non-radiative escape channels and relaxation 
processes need to be taken into account. For thicker GaTe films, much stronger PL than predicted 
by the model is observed, which probably signifies suppression of additional non-radiative escape 
that free excitons experience in relatively thin films. It is also notable that in GaTe the PL lifetime 
for localized states is shorter than in GaSe, which may be due to non-radiative escape. Such 
processes are not included in the model. Another reason could be deviation from the effectively 
’ballistic’ transport that we assume leads to the carrier escape into regions (1), and its replacement 
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for the larger thicknesses with a slower ’diffusion’ process leading to slower non-radiative escape. 

Assuming that the measured PL rise-times of ~ 20 ps are close to T re u we can estimate char¬ 
acteristic non-radiative times x nr \ and z nr 2 for several limiting cases. For example for GaSe we 
obtain the following values using T, e /Wi=1000 nm and t, y ./M2=K)0 nm: for hfu m = 10 nm t,„-i=0.05 
ps, for hfu m =30 nm t,„ i=0.15 ps, for h/u rn = 100 nm x nr 2 =7 ps. For both h fn m of 10 and 30 nm, the 
non-radiative decay occurs on a sub-picosecond time-scale. Here, the non-radiative escape is by a 
factor of the order of 100 faster than relaxation into the light emitting states, which is consistent 
with a low quantum yield of 10 3 reported previously for thin films of M 0 S 2. 9 For the middle 
region of a 100 nm film, the characteristic non-radiative escape time is comparable with the relax¬ 
ation time into the states giving rise to PL, thus a much higher quantum yield can be expected for 
films of this thickness. 

We also note that in the case of region (2), 112 should be of the order of a typical thermal velocity 
of an exciton, v t i v By using T = 10 K and the exciton mass of 0.1 (or 0.2) 39 of the free electron 
mass we obtain T re /v f /, ~ 1200 (or 900) nm for T re i ~20 ps, very similar to the order of magnitude 
predicted by the model for the product x re iU 2 . 

Conclusions 

We have investigated thin films of GaTe and GaSe prepared by mechanical exfoliation from bulk 
crystals, deposited on Si02 substrates and capped with a thin layer of The study of optical 

properties of the thin films has been conducted by means of low-temperature cw and time-resolved 
micro-PL techniques. The most pronounced property observed is a significant reduction of cw PL 
intensity by up to 10 5 for thin films of thicknesses about 10 nm compared with the films of 150-200 
nm. No measurable PL was observed in films thinner than 7 nm. Except for PL decrease, no other 
clear trends, including size-quantization effects and PL life-time modifications as a function of the 
film thickness were observed in cw and time-resolved PL. We argue that the reduction of quantum 
yield occurs due to the non-radiative processes associated with the surface states. The theoretical 
model that we develop differentiates between the fast non-radiative carrier escape to surface traps 
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in the thin 15-20 nm layers near the film surface, and slower decay in the middle region of the 
film, where the carriers first move to the surface layers and then decay. The model accounts for the 
change in the character of the PL decay for thin films with thicknesses less than 30-40 nm, where 
we expect only the fast direct relaxation to surface traps. 

We do not refer in our model and interpretation to the direct-to-indirect bandgap transition 
found for thin InSe films in Ref.— The observed strong PL decrease is not accompanied with any 
significant variation of the PL lifetime, which would most likely accompany such a change in the 
band structure. This implies that the PL reduction in our case is not associated with this type of 
transition, and instead the reduction in the radiative recombination contribution is likely due to an 
increase in the non-radiative decay rate in thin films. 

The strong non-radiative decay processes occur in the studied films despite complete encapsu¬ 
lation in Si 3 N 4 and Si 02 providing partial protection of the surface from chemical interactions with 
the ambient atmosphere. This emphasizes the need for development of novel surface passivation 
for III-VI films, possibly involving oxygen free substrates and encapsulation in additional layered 
materials such as boron nitride. Using this fabrication methods, the large family of III-VI materials 
may show weaker dependence of quantum yield on film thickness, which will allow them to play 
important role as building blocks in van der Waals heterostructures. It is likely that in the near fu¬ 
ture fabrication of such heterostructures will be carried out in oxygen- and water-free atmospheres 
and will include carefully designed passivating layers. 
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Figure 1: (a) AFM image of a GaSe thin film having a single unit cell thickness. This film is 
obtained by mechanical exfoliation, (b) Cross-sectional plot along the black line in (a). 
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Figure 2: Low-temperature cw PL results for GaTe thin films, (a) PL spectra for films of various 
thicknesses. hfu, n is marked on the plots. On the plot for the film with hfn m = 49 nm FE marks 
the free exciton peak observed at similar energies for all films on the figure except the one with 
hfiim =8 nm. (b) Symbols show experimentally measured integrated PL (for 7=1 OK) for the cw 
laser excitation power of 2 mW at 532 nm (2.33 eV). Dashed line shows expected variation of 
PL following the change in the absorption of the thin film assuming constant quantum efficiency. 
Solid line shows the results of calculations using the model discussed in text. 
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Figure 3: Low-temperature cw PL results for GaSe thin films, (a) PL spectra for films of various 
thicknesses, (b) Symbols show experimentally measured integrated PL (for 7=1 OK) for the cw 
laser excitation power of 2 mW at 532 nm (2.33 eV). Dashed line shows expected variation of 
PL following the change in the absorption of the thin film assuming constant quantum efficiency. 
Solid line shows the results of calculations using the model discussed in text. 
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Figure 4: (a) Observation of the Stokes shift of localized exciton states in cw PL of a 28 nm thick 
GaSe films. Excitation with a cw laser at 532 nm (2.33 eV) is used at T=10K. (a) Observation 
of the free exciton feature in PL of a 42 nm thick GaSe films. The time-integrated PL is shown 
measured for pulsed laser excitation at 420 nm (2.95 eV) with a power of 1 mW. A free exciton 
peak having a linewidth of ^25 meV is observed at 2.06 eV. 
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Figure 5: Time-resolved PL data for GaTe (a,b) and GaSe (c,d) thin films measured at T=10K. 
Note different time-scales for the GaTe and GaSe data, (a) and (c) show streak-camera scans, 
(b) and (d) show PL traces measured at 1.76 (gray) and 1.71 eV (red) for GaTe and 2.055 (gray) 
and 2.035 eV (red) for GaSe. Blue lines in (b) show fitting with a single exponential functions as 
described in text. 
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Figure 6: Diagrams illustrating the theoretical model for the dependence of the cw PL intensity on 
the film thickness, (a) Processes included in the model (see text and supplementary information for 
more details): optical excitation of e-h pairs; relaxation into the non-radiative traps with a time T nr \ 
or T nr 2 depending on the thickness of the film; relaxation with a time T re i into the light-emitting 
states denoted as ’PL states’; PL emission from these states with a time t pl- (b) and (c) shows two 
types of thin films with thicknesses above and below the critical thickness of 21t.q. In (a) we show 
a relatively thick film with a thickness of hju m > 2Hq. This film has three regions: two regions of 
thickness ho near the film surfaces where fast carrier relaxation to surface traps occurs (shown in 
pink), and a ’normal’ region of thickness hfn m — 2ho in the middle of the film (blue) where non- 
radiative processes are weaker and occur through carrier escape into the surface regions. In (b) we 
show a thin film with hfu m < 2Iiq, where a region of one type only exists, where fast non-radiative 
carrier decay occurs. 
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Optical properties of two-dimensional gallium chalcogenide 
films: Supplementary Information 


Abstract 

Here we present additional details on the rate-equation model describing the thickness depen¬ 
dence of photoluminescence intensity in GaSe and GaTe thin films. 


Modeling 

In order to describe the observed trend of photoluminescence (PL) intensity as a function of film 
thickness hju, n in GaSe and GaTe thin films we have developed a simplified rate equation model 
described below. It describes the behavior of photo-excited carriers in terms of populations of e-h 
pairs, i.e. electrons and holes are not treated separately for simplicity. 

As shown in Fig fSTt a). in the model we consider the following processes leading to light emis¬ 
sion: optical excitation of e-h pairs (uniformly across the full thickness of the film); relaxation into 
the non-radiative traps with a time T nr \ or z nr 2 depending on the thickness of the film, as explained 
below; relaxation with a time t re i into the light-emitting states denoted in FigIS 1 la) as ’PL states’; 
PL emission from these states with a time Xpp. The absence of a clear dependence of the PL decay 
times on the film thickness allows us to use a thickness-independent constant T pp and neglect any 
other decay mechanisms for the e-h population of the light-emitting states. We also assume that 
the mechanism of relaxation into the ’PL states’ is not sensitive to the film thickness and can be 
described with a thickness independent time x re i. 

As shown in Fig lSTT b) and (c) we assume that the film is divided into three regions: (1) two 
regions of thickness ho near the film surfaces where fast carrier relaxation to surface traps occurs 
leading to non-radiative carrier escape; (2) a ’normal’ region of thickness h fu m — 2ho in the central 
part of the film where carriers do not decay non-radiatively, but can escape into the regions (1), 
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Figure SI: Diagrams illustrating the theoretical model for the dependence of the cw PL intensity 
on the film thickness, (a) Processes included in the model (see text and supplementary information 
for more details): optical excitation of e-h pairs; relaxation into the non-radiative traps with a time 
T nr \ or t ,„-2 depending on the thickness of the film; relaxation with a time x re i into the light-emitting 
states denoted as ’PL states’; PL emission from these states with a time T>/,. (b) and (c) shows two 
types of thin films with thicknesses above and below the critical thickness of 2ho. In (a) we show 
a relatively thick film with a thickness of hjn m > 2/io . This film has three regions: two regions of 
thickness ho near the film surfaces where fast carrier relaxation to surface traps occurs (shown in 
pink), and a ’normal’ region of thickness hfu, n — 2/?o in the middle of the film (blue) where non- 
radiative processes are weaker and occur through carrier escape into the surface regions. In (b) we 
show a thin film with hfu m < 2ho, where a region of one type only exists, where fast non-radiative 
carrier decay occurs. 


where they undergo non-radiative decay. Here the ho value may be associated with a depletion 
depth or an average surface trap radius. S1 It is assumed to be independent of the film thickness and 
may be different for films made of different materials. 

In region (1), in first approximation the average time it takes for the carrier/e-h pair to escape 
non-radiatively is proportional to half the thickness of region (1) (can be understood as the average 
time for the carrier to reach the surface or as the overlap of the wavefunction of the carrier and 
the trap). This half-thickness is given by ho/2 for the two regions adjacent to the surface for 
film thicknesses hfu, n > 2ho, and is hfu m /4 for hfu m < 2ho, when the two surface regions merge. 
Thus we introduce a non-raditive decay time in the regions of type (1) as T nr \ — {ho/2)/u\ for 
h film > 2/?o, and x nr \ — (/?/,/,„/4) /u\ for hfu, n < 21iq. Here u\ is a constant with the dimensions of 
meter per second. 

In region (2), the non-radiative escape time reflects the average time it takes for a carrier or 
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an e-h pair to reach any of the regions (1). The underlying mechanism for this process may be 
depletion and band-bending expected at the film surface leading to charge separation and non- 
radiative decay. S1 We assume that once the carrier or e-h pair has reached region (1) it escapes 
non-radiatively. In first approximation, the average time it takes a carrier/e-h pair to reach region 
(1) is proportional to half the thickness of region (2), h fu m /2 — ho. Thus we introduce a non-raditive 
decay time in region (2) as t ,„2 = [hfu m /2 — ho)/U 2 - Here U 2 is a constant with the dimensions of 
meter per second, effectively corresponding to the average carrier velocity in region (2). Notably, 
in films with hfu m < 2ho, region (2) does not exist, and the only non-radiative escape that we 
consider is due to the escape to traps as described above for region (1). 

In equations below, we denote e-h pair populations in regions (1) and (2) as N\ and Afe, respec¬ 
tively, and the population of the ’PL states’ as Ni oc . Appropriately, A3 — 0 for hfu m < 21iq. In both 
regions, e-h pairs relax with the time t re i into the ’PL states’ giving rise to radiative recombination. 
Non-radiative escape from the ’PL states’ is neglected. This assumption is based on the fact that 
in most cases PL is observed from the tightly-bound localized states. Ni oc states decay radiatively 
with a characteristic time T PP , so that the PL intensity I PP = N[ oc /% pp . Saturation of ’PL states’ 
is neglected as the data reported in Figs JS2IS~3l is measured for the laser powers where saturation 
effects are unimportant. The light absorption coefficient, a is assumed to be the same in both types 
of regions (1) and (2). For the small thicknesses of the films that we consider, the e-h pair genera¬ 
tion rate is equal to Pa multiplied by the thickness of the regions (1) or (2) as given below. Here 
P is the e-h pair generation rate due to the laser excitation rescaled to account for the PL detection 
efficiency. 

The following rate equations can thus be introduced for the case hjn m > 2ho'. 

dN\ / dt — 2hoPa — Ni jt re i — N\ / x nr \ 

JAA/ dt — (hfii m — 2ho)Pa —N 2 /t re i —N 2/% tr 2 

dNioc/dt N\ / T re i -\- N 2 / T r et N[ oc /x PP 

For the case hfu m < 2/iq we will have: 
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dN\/dt P hfilm CC N\/t re i Ni/x nr i 
N 2 =0 

dN\\ oc j dt N\ / T re i N,oc/tpL 

In the steady-state case we obtain the following expression for PL intensity for hfu m > 2h 0 : 


IpL — 2h()PCC/ (Id - 'tyel I %irl ) {hfilm 2Ilo)Paj (Id - X re ij T nr2)i (1) 


where T nr \ = T|/? 0 /2 and x nr2 = x 2 {hfn m /2 - /to), and for h f u m < 2h 0 


IpL hfil m POC / ( 1 “I - Xygl I 


( 2 ) 


where T m -, = {h film /4)/u h 

By substituting these expressions for the non-radiative decay times directly we obtain for PL 
intensity for hfu, n > 2/iq: 


Ipl = 2 h 0 Pa/{ho + 2t re ju\) + {hfum — 2ho) 2 Pa/{hfu m — 2ho + 2x re iU2). (3) 

For hfum < 2/7o we get: 

Ipl = hf Um P a / (hfum + 4 T re / u i), (4) 

Results of the modeling 

Results for GaSe thin films 

According to Eqs. ??, ??, we have four parameters, with which we describe the observed depen¬ 
dence of Ipl on the film thickness: ho, the products Pa, x,ei u \ and V/ M 2 - 

Further analysis of Eqs. ??,?? shows that these parameters are not completely independent. 
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Figure S2: Low-temperature cw PL results for GaSe thin films. Symbols show experimentally 
measured integrated PL (for T =10K) for the laser excitation power of 2 mW (laser excitation at 532 
nm is used). Dashed line shows expected variation of PL following the change in the absorption of 
the thin film. Solid line shows the results of calculations using the model discussed in text. 

This becomes obvious if we note that for small hfu, n where we expect fast non-radiative process 
(i.e. large x re iu\) we get h y,/ m -C i , and Eq.?? reduces to a simple parabolic dependence: 
I PL h 2 r ilm Pd/ (4z re iu\). This function is independent of individual parameters Pa and z re iu\, and 
depends on their ratio only. 

For large hfu m where we expect slower non-radiative process we get ( hfu m — 2ho) of the order 
of 2x re iU2, and a simple approximation can be used for Eq. ?? giving Ipi ~ (hfu m — 2ho)Pa. From 
here we find that Pa should be of the order of 10 5 nm 1 s 1 . 

The two functions described by Eqs. ??, ?? and especially their asymptotics have significantly 
different form, effectively with different ’slopes’ as seen e.g. in FigJS2]for the case of GaSe. For 
GaSe data in FigJS2l it is relatively easy to define the position of a characteristic kink around 
hfn m =30 nm, where the behavior changes. We use this observation to fix the value of ho to 15 nm. 

In the next step, we have explored a range of fitting functions providing a reasonable agreement 
with the experiment. We conclude that our model can reliably predict the order of magnitude of 
the fitting parameters, but not the actual value. So we find that z re iu\ is of the order of a 1000 nm 
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Figure S3: Low-temperature cw PL results for GaTe thin films. Symbols show experimentally 
measured integrated PL (for T =10K) for the laser excitation power of 2 mW (laser excitation at 532 
nm is used). Dashed line shows expected variation of PL following the change in the absorption of 
the thin film. Solid line shows the results of calculations using the model discussed in text. 

and T re iU 2 is of the order of a 100 nm. 

As an illustration for the GaSe data, in FigJS2]we show a fitting with Pa = 3.5 x 10 5 nm ! s 
T re [U\ = 2900 nm and t ,- c /«2 = 150 nm. The figure also shows the variation of PL intensity as 
would be expected from the change in absorption only (dashed line). This would correspond to a 
thickness independent non-radiative processes occurring with the same rate in the whole volume of 
the film. This curve is described by the expression Ipi = IcaSe[\ ~ ^xp{~ a GaSehfiim)\, where IcaSe 
is the PL intensity for films with hfu m ~200 nm and the absorption coefficient a= 1000 cm .— 

Results for GaTe thin films 

We used the same fitting procedure for GaTe with the results shown in FiglS3l Although the ’kink’ 
in the thickness-dependence of Ipp for GaTe is not as pronounced as for GaSe, it is still observable 
at a slightly higher value of 2fio=40 nm. The model describes adequately the data for hfu m up to 
around 90 nm, but fails to reproduce the very steep enhancement of PL for thicker films. This 
probably occurs due to the contribution of the free excitons to the PL signal. It is probable that 
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there is an additional strong non-radiative escape mechanism for free excitons in thin films, which 
is overcome in thicker films with the result that the PL intensity for thick films grows faster than the 
model predicts. Another reason could be deviation from the effectively ’ballistic’ transport that we 
assume leads to the carrier escape into regions (1), and its replacement for the larger thicknesses 
with a slower ’diffusion’ process. 

In a similar way to GaSe, we confine ourselves to determination of the order of magnitude 
only for the main parameters, which are similar to GaSe: we find that T re iU\ is of the order of a 
1000 nm and x re iU 2 is of the order of a 100 nm. As an illustration for the GaTe data, in FiglS3lwe 
show a fitting with Pa = 5.0 x 10 6 nm^'s' 1 , x re i/X\ = 3800 nm and T re //T 2 = 200 nm. The figure 
also shows the variation of PL intensity as would be expected from the change in absorption only 
(dashed line). This curve is described by the expression I PL = Ic a Te[ 1 — exp(—<XGaTehfu m )], where 
iGaTe is the PL intensity for films with h fn m ~200 nm. The absorption coefficient oy ; - a 7>=5000 
cm -1 is used according to Ref. S2 
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